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Abstract 

This paper presents measurements from the ATLAS experiment of the forward-backward 
asymmetry in the reaction pp —» Z/y* —> 1 + 1~, with I being electrons or muons, and 
the extraction of the effective weak mixing angle. The results are based on the full set 
of data collected in 2011 in pp collisions at the LHC at sfs = 7 TeV, corresponding 
to an integrated luminosity of 4.8 fb -1 . The measured asymmetry values are found to 
be in agreement with the corresponding Standard Model predictions. The combination 
of the muon and electron channels yields a value of the effective weak mixing angle of 
sin 2 = 0.2308 ± 0.0005(stat.) + 0.0006(syst.) ± 0.0009(PDF), where the first uncertainty 
corresponds to data statistics, the second to systematic effects and the third to knowledge 
of the parton density functions. This result agrees with the current world average from the 
Particle Data Group fit. 
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1 Introduction 


The vector and axial-vector couplings in the neutral current annihilation process qq —* Z/y* —> t + ('~ 
lead to a forward-backward asymmetry Afb in the polar angle distribution of the final state lepton t~ with 
respect to the quark direction in the rest frame of the dilepton system. This paper presents measurements 
of the forward-backward asymmetry in electron and muon pairs from Z/y* boson decays and the extrac¬ 
tion of the weak mixing angle by the ATLAS experiment. The results are based on the full set of pp 
collision data collected in 2011 at the LHC at a centre-of-mass energy of s[s = 7 TeV, corresponding to 
an integrated luminosity of 4.8 fb -1 . 

The differential cross section for the annihilation process can be written at leading order as 


dcr 

dfcos 9) 


4na~ 


-A( 1 + cos 2 9) + B cos 9 


( 1 ) 


where a is the fine-structure constant, Vs is the centre-of-mass energy of the quark and anti-quark, and 
9 is the angle between the lepton and the quark in the rest frame of the dilepton system. The coefficients 
A and B arc functions of VS and of the electroweak vector and axial-vector couplings. In the case that 
the dilepton system has non-vanishing transverse momentum, p t, the four-momentum of the incoming 
(anti-)quark is not known, as it is no longer collinear with the incoming beams. The impact of this effect 
on the asymmetry measurement is minimized by choosing a particular rest frame of the dilepton system, 
the Collins-Soper (CS) frame [1], in which the angle between the lepton and the quark, 9* cs , is calculated. 
The sign of cos 9‘ cs is defined with respect to the direction of the quark, which is, however, ambiguous in 
pp collisions. It is therefore chosen by measuring the longitudinal boost of the final-state dilepton system 
in the laboratory frame, and assuming that this is in the same direction as that of the quark in the initial 
state. This assumption leads to a fraction of events with wrongly assigned quark direction, which causes a 
dilution of the observed asymmetry. The probability of correct quark direction assignment increases with 
the boost of the dilepton system, thus reducing the dilution for dileptons produced at large rapidities. With 
this assumption, cos 9*, s can be written as a function of the lepton momenta in the laboratory frame, 


Pz,tC 2 (PiP 2 ~PiPV 
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( 2 ) 


with 


Pi = ~p( E i ± Pz,i), 

V2 

where E is the energy and p z the longitudinal momentum of the lepton (i = 1) and anti-lepton (i = 2). 
The variables p,j(, nice, and pjj/ denote the longitudinal momentum, invariant mass and transverse 
momentum of the dilepton system, respectively. The first factor in eq. 2 defines the sign of cos 9* cs 
according to the longitudinal direction of flight of the dilepton system, as discussed above. The events 
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with cos (9*- s > 0 are classified as forward (F), while those having cos 6* cs < 0 are classified as backward 
(B). The asymmetry Afb is then defined as 


Afb 


o~f - o~ B 
CTp + o- B ’ 


(3) 


where c rp and <jb are the cross sections for the respective forward and backward configurations. At 
leading order, the second term in eq. 1, Bcos 9, describes the asymmetry Afb- This analysis measures 
Afb as a function of the invariant mass of the dilepton system. The results, which are presented in 
section 5, include the detector-level values, as well as the corrections needed to take into account detector 
effects and dilution. 

Several Standard Model parameters can be extracted from the dependence of the Afb values on the invari¬ 
ant dilepton mass. One of these is the electroweak mixing angle sin 2 Ow, which is defined at tree level as 
1 - Depending on the renormalisation scheme, higher-order loop corrections may modify this 

relation. This analysis extracts the effective leptonic weak mixing angle sin 2 0^ [2, 3] at the my scale 
from the detector-level Afb values. 

The effective weak mixing angle is related to the electroweak vector coupling g v via 

9v = VPf{ T j ~ 2 Qf sin2 #e//) , with sin 2 8 eff = Kf sin 2 6 W , 

where the electroweak radiative corrections to the tree-level couplings are absorbed into the fermion- 
dependent factors Kf and py, Tl is the third component of the weak isospin and Qf the electric charge of 
the fermion /. Using this definition of the effective weak mixing angle, the coupling retains its tree-level 
form multiplied by the additional factor yfpj. The relationship between the leptonic and quark sin 2 $ e yy 
can be approximated as a flavour-dependent shift in the leptonic sin 2 0 e yy [3]. Although the sin 2 0 e yy value 
from 6-quarks differs the most from the one from leptons, only a few percent of the events in this analysis 
come from initial-state 6-quarks. In particular, the effect of the quark sin 2 8 c fj on the measured Afb is an 
order of magnitude smaller than the effect of the leptonic sin 2 0 e //- This analysis therefore measures the 
leptonic sin - 0 e yy, denoted by sin - 8 c j j in the following. Its value is extracted from the measured Afb as 
a function of the invariant mass of the dilepton system by comparing it to MC predictions produced with 
varying values of the weak mixing angle. Details are given in section 6. 

The most precise measurement of sin 2 8^ comes from the combination of results from the LEP and SLD 

experiments [3]. Those studies yield an average leptonic sin 2 8^ = 0.23153 ± 0.00016. The two most 
precise single measurements are extracted from the forward-backward asymmetry in 6-quark final states, 
Ap*, at LEP (sin 2 8^ = 0.23221 ± 0.00029) and from the leptonic left-right polarization asymmetry, 

Alr, at SLD (sin 2 6^j = 0.23098 ± 0.00026). These two values differ by approximately three standard 
deviations. More recently, the CDF [4] and DO [5] experiments at the Tevatron and the CMS [6, 7] 
experiment at the LHC have also measured sin 2 d^ff- The CDF (DO) measurement was performed using 
Z —* pp (Z —* ee) events from pp collisions, and the CMS measurements were performed using Z —> pp 
events from pp collisions. These results are compared to those from this analysis in section 6. 

The value of Afb at the peak of the Zjy* resonance {mu — my), A |B , can be written as a function of the 
asymmetry parameters Ac and A q , 




( 4 ) 
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with £ (q) denoting the leptons (quarks) in the final (initial) state. The parameters A/ and A q are directly 
related to the electroweak vector and axial-vector couplings, as described in section 6.3. The most precise 
measurements of the electron and muon asymmetry parameters were performed by SLD [3], yielding 
A e = 0.15138 ± 0.00216 and A p = 0.142 ± 0.015. The precision of the A q measurement is dominated 
by the statistical uncertainty, thus making it an interesting parameter to measure with the large number of 
Z —> ini events produced at the LHC. The A fl result from this analysis is presented in section 6.3. The 
determination of A fJ in the LEP/SLD results is entirely based on asymmetry measurements in electron 
and muon final states without any assumptions on the involved Af. In contrast, the determination of A fl 
presented here uses the Standard Model prediction of A q . 


2 The ATLAS detector 

The ATLAS detector [ 8 ] is a general-purpose detector installed at the LHC [9] at CERN. The detector 
subsystem closest to the interaction point, the inner detector (ID), provides precise position and mo¬ 
mentum measurements of charged particle trajectories. It covers the pseudorapidity range 1 \q\ < 2.5 and 
provides full azimuthal coverage. The ID consists of three subdetectors arranged in a coaxial geometry 
around the beam axis: the silicon pixel detector, the semiconductor microstrip detector and the straw- 
tube transition-radiation tracker. A solenoidal magnet generates a 2 T magnetic field in which the ID is 
immersed. 

Electromagnetic calorimetry in the region |//| < 3.2 is based on a high-granularity, lead/liquid-argon (LAr) 
sampling technology. Hadronic calorimetry uses a scintillator-tile/steel detector in the region |//| <1.7 
and a copper/LAr detector in the region 1.5 < |//| < 3.2. The most forward region of the detector 
(3.1 < I/ 7 I < 4.9) is equipped with a forward calorimeter, measuring both the electromagnetic and hadronic 
energies using copper/LAr and tungsten/LAr modules. 

A large muon spectrometer (MS) constitutes the outermost part of the detector. It consists of three large 
air-core superconducting toroidal magnet systems (each with eight coils): one barrel providing a field of 
about 0.5 T and two endcaps each providing a field of about 1 T. The deflection of the muon trajectories 
in the magnetic field is measured in three layers of precision drift tube chambers for \q\ < 2. In higher 
I77I regions (2.0 < |//| < 2.7), two layers of drift tube chambers are used in combination with one layer of 
cathode strip chambers in the innermost endcap wheels of the MS. Three layers of resistive plate chambers 
in the barrel (|//| < 1.05) and three layers of thin gap chambers in the endcaps (1.05 < |?/| < 2.4) provide 
the muon trigger and also measure the muon trajectory in the non-bending plane. 

A three-level trigger system is used to select events in real time. A hardware-based Level-1 trigger uses 
a subset of detector information to reduce the event rate to a design value of at most 75 kHz. The rate 
of accepted events is then reduced to about 300 Hz by two software-based trigger levels, Level-2 and the 
Event Filter. 


1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points 
upward. Cylindrical coordinates are used in the transverse plane, <t> being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as rj = - In tan(0/2). 
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3 Signal and background modelling 


Monte Carlo (MC) simulated event samples used to model signal and background processes are generated 
and passed through the ATLAS detector simulation [10], based on the GEANT4 toolkit [11]. Simulated 
events acquire weights such that the resulting distributions match the ones observed in the data for the 
following variables: the average number of interactions per bunch crossing, the z coordinate of the in¬ 
teraction vertex, the lepton energy/momentum scale and resolution, and the trigger, identification and 
reconstruction efficiencies. 

The Z/y* production is detected by the emission of charged lepton pairs, ee or pp. The contribution from 
Z/y* —> rr followed by r decays to electrons or muons is considered as background and subtracted from 
the signal. Signal samples are generated with PYTHIAv6.4 [12] with the MSTW2008LO [13] parton 
distribution functions (PDFs) and a value of sin 2 9^j?* = 0.232 for the effective weak mixing angle. Final- 
state radiation from QED is taken into account using PHOTOS [14] in the exponentiated mode with multi¬ 
photon emission. For the sin 2 6^ measurement, where sensitivity to PDFs is expected to be significant, 
additional PDF sets are also used, including one specifically prepared for this analysis, based on the 
ATLAS-epWZ12 PDFs [15]. Details are given in section 6. The cross section is calculated at next-to-next- 
to-leading order (NNLO) in the strong coupling using PHOZPR [16] with MSTW2008 NNLO PDFs. The 
ratio of this cross section to the leading-order (LO) cross section is the m^-dependent /("-factor, applied 
to the generated signal for all plots shown in the following. However, the main observable described here 
(Arb) is not affected by this LO-to-NNLO rescaling. The impact of higher-order corrections in and a em 
on Afb and on sin 2 is assessed using the HORACEv3.1 [17], MCFMv6.6 [18] and POWHEGvl [19] 
generators, as described in section 5. The POWHEG simulation is combined with PYTHIA 6.4 for 
showering and hadronization. 

Background contributions containing prompt isolated electron or muon pairs are estimated using Monte 
Carlo simulation normalized using the best available cross section prediction at (N)NLO. The background 
from Z/y* —> rr is also generated using PYTHIAv6.4. Diboson ( WW, WZ, and ZZ) samples are gen¬ 
erated with HERWIGv6.510 [20, 21] using the MRSTMCal PDFs [22], Pair-production of top quarks is 
generated with MC@NLOv4.01 [23, 24] using the CTQE6L1 PDFs [25], combined with HERWIG for 
showering and hadronization. 

The contributions from multi-jet and VT+jcts background events containing non-isolated leptons from 
heavy-flavour decays and hadrons misidentified as leptons are estimated using data-driven techniques, 
as described in section 4.3. Since the contribution from VV'+jcts is found to be a small fraction of the 
multi-jet background over the whole invariant mass range, the term ‘multi-jet background’ is used in the 
following to denote the sum of these contributions. 


4 Event reconstruction and selection 

The analysis uses pp data collected in 2011, corresponding to an integrated luminosity of 4.8fb _1 for 
the electron channels and 4.6fb _1 for the muon channel. All events analysed were acquired under good 
operating conditions of the ATLAS detector. Events in the electron channels passed the single electron 
trigger, with an electron Ej > 20 or 22 GeV (depending on the instantaneous luminosity). Events in the 
muon channel passed the single muon trigger with a muon pj threshold of 18 GeV. The presence of a 
reconstructed collision vertex with at least three tracks with pj > 400 MeV is required. For the muon 
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channel, there is an additional requirement that the longitudinal position of this vertex be within 200 mm 
of the nominal interaction point. 


4.1 Electron reconstruction 

This analysis uses electrons in two distinct regions of the detector: the central region (|//| < 2.47) where 
tracking information is available, and the forward region (2.5 < \rj\ < 4.9), where the electron reconstruc¬ 
tion relies only on information from the calorimeter. The inclusion of electrons in the forward region 
allows the reconstruction of events where the Z/y* candidates are emitted at large rapidity, thus reducing 
the effect of dilution due to the imperfect knowledge of the direction of the initial state quark. 

For both the central and forward electrons, the reconstruction begins with identifying energy deposits in 
the calorimeters consistent with electromagnetic showers. Electron candidates in the central region arc 
matched to a track reconstructed in the ID. A transverse energy requirement, Ej > 25 GeV, is applied to 
both the central and forward candidates. Electron candidates in transition regions between the barrel and 
endcap calorimeters (1.37 < \r/\ < 1.52) and between the endcap and forward calorimeters (3.16 < |//| < 
3.35) are excluded from this analysis. 

The central candidates must satisfy either ‘medium’ or ‘tight’ identification criteria, based on shower 
shape and track quality variables [26] optimized for the 2011 data [27]. Forward electron candidates must 
satisfy similar medium quality criteria optimized specifically for forward electrons. When combining one 
central and one forward electron there is an additional requirement that the central electron be isolated. It 
is implemented by requiring that the transverse energy deposition in the calorimeter within a cone defined 

by AT? = ^J{A(p) 2 + (A?/) 2 = 0.2 around the electron candidate be less than 5 GeV, excluding the electron 
candidate itself. The reconstruction efficiencies in simulated events are corrected to match the measured 
efficiencies [27], Selected events consist of either two medium candidates in the central region (central- 
central, referred to as CC) or one tight central electron candidate and one medium forward electron 
candidate (central-forward, referred to as CF). In the CC electron channel, the electrons are required to 
have opposite charges. This requirement is not applicable in the CF electron channel, since the forward 
electron has no charge information. The effect of charge misidentification is found to be negligible in 
both the CC and CF electron channels [27], 


4.2 Muon reconstruction 

Muons identified as ‘combined muons’ by the reconstruction and identification algorithms [28, 29] are 
used in this analysis. They arc reconstructed by associating and combining two independently reconstruc¬ 
ted tracks, one in the ID and one in the MS. Combined muons are required to have transverse momentum 
pj > 20 GeV, and must lie within |t/| < 2.4. The ID tracks associated with the muons must satisfy quality 
requirements on the number of hits recorded by each subdetector [28], To reject muons from cosmic rays, 
the longitudinal coordinate of the point of closest approach of the track to the beamline is required to be 
within 10 mm of the collision vertex (see section 4.3). Rejection of multi-jet background is improved by 
requiring the muons to be isolated. The isolation parameter is the relative momentum isolation, defined 
as the sum of the p \ of all other tracks within a cone of AT? = 0.2 around the muon track, divided by 
the muon p-\\ Y Pj"‘ k I Pj < 0.1. The kinematic variables of the muons are measured by the ID, in order 
to minimize the impact of residual misalignments between the ID and the MS. This choice also reduces 
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the impact of muon bremsstrahlung in the calorimeter on the measurement. Charge misidentification for 
muons is very low, with negligible effect on this analysis. 


4.3 Event selection and background estimation 

Events must contain two oppositely charged leptons in the muon and CC electron channels or one central 
electron and one forward electron in the CF electron channel. In the muon and CC electron channels, 
dilepton pairs with invariant masses up to 1000 GeV are used. In the CF electron channel, the Afb meas¬ 
urement is performed only for dilepton masses up to 250 GeV, because the background dominates at larger 
masses, leading to sizeable systematic uncertainties. 

Contributions of different background sources are estimated using either simulation or data-driven tech¬ 
niques. For dibosons (ZZ, WZ, WW), Z/y* —> rr and tt, contributions are estimated using simulation. The 
dominant background for the muon and CC electron channels, across the whole invariant mass range, is 
that due to tt events. Background from Z/y* —» tt production (followed by r —> tv) populates the low end 
of the dilepton invariant mass distribution. For the electron channels, the multi-jet background is estim¬ 
ated using a combination of data-driven techniques. In the CC electron channel, the reverse identification 
method [30] is used for dilepton invariant masses below 125 GeV, while the fake-factor method [31], is 
employed for higher invariant masses. An overlap region is defined between 110 and 200 GeV, where 
the estimates from both methods are compared and a scale factor for the reverse identification estimate is 
determined using the fake-factor result. Since the CF electron channel only extends to a dilepton invari¬ 
ant mass of 250 GeV, only the reverse identification method is used. In the muon channel, the multi-jet 
background is estimated from data in a control region defined by inverting the muon isolation cut. The 
numbers of events in the control and signal regions observed in MC simulation are then used to transfer 
the data distribution from the control region to the signal region. 

Figures 1 and 2 show the met and cos (9* ;s distributions of events in the three channels. The total numbers 
of selected events arc 1.2xl0 6 , 0.35xl0 6 and 1.7xl0 6 for the CC electron, CF electron and muon channels 
respectively. In the region close to the Z peak, the background contamination is estimated to be less than 
1% for the muon and CC electron channels, and about 5% for the CF electron channel. The background 
contributions in the muon and CC electron channels increase to about 5% and 16% in the low- and high- 
mass regions, respectively. The CF electron channel has a background contamination of about 30% in 
the low-mass region. Agreement between data and simulation is observed within the uncertainties over 
the whole invariant mass range and also in the cos 6* cs distributions. These uncertainties contain both 
the statistical and systematic components and include the effects of multiple pp collisions occurring in 
the same or in neighbouring bunch crossing (pileup), energy/momentum scale and resolution, trigger 
efficiency, misalignment of the inner detectors, data-driven background estimates, and PDFs. Details are 
given in section 5. 

Figure 2(c) highlights the cos 9* cs distribution for the CF electron channel to better illustrate the reduced 
impact of dilution: the forward-backward asymmetry is large enough to be observed directly from the plot. 
Some differences between data and simulation are observed in the lowest and highest bins in cos 0* s . As 
a cross check, the analysis was repeated excluding the bins in question and the impact on the Afb and 
sin 2 (T^results was found to be negligible. 




Figure 1: Dilepton invariant mass distributions obtained from the event selections described in the text, for the (a) 
CC electron, (b) CF electron and (c) muon channels. Data are shown by open circles and the total expectation is 
shown as a line with a band representing the total uncertainty (statistical and systematic added in quadrature). The 
data-driven estimate for the multi-jet background and the simulation-based estimates for all other backgrounds are 
shown by the shaded areas. 
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(c) (d) 

Figure 2: Distributions of the cosine of the polar angle in the Collins-Soper frame (cos f)’ cs ) obtained from the event 
selections described in the text, for the (a) CC electron and (b) muon channels. The corresponding distribution for 
the CF electron channel is shown using both (c) a linear and (d) a logarithmic scale. Data are shown by open circles 
and the total expectation is shown as a line with a band representing the total uncertainty (statistical and systematic 
added in quadrature). The data-driven estimate for the multi-jet background and the simulation-based estimates for 
all other backgrounds are shown by the shaded areas. 
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5 Measurement of A 


FB 

For each invariant mass bin, the Afb value is obtained from the corresponding cos d*, s distribution by 
measuring the numbers of forward and backward events: 

N cos WW>() — N cos (f rK <0 

A fb =--—. 

^cos6J s >0 "I" A^ CO s0* s <O 

For comparison, expected Afb values are calculated from both the PYTF1IA and POWHEG samples de¬ 
scribed in section 3. Background contributions are subtracted from the number of forward and backward 
events measured at detector-level. Some background contributions, such as multi-jet events, display no 
asymmetry and hence dilute the measured asymmetry. Other background contributions, such as tt, display 
an asymmetry. The detector-level asymmetry values after background subtraction (Aj?g as ) in the electron 
and muon channels are shown in figure 3 as a function of the invariant mass of the lepton pair 2 . Good 
agreement between data and simulation is observed. Figure 4 shows the same information in a narrower 
mass range around the Z pole. 

The asymmetry values A j™ as are unfolded from detector level to particle level (A°g S ), to allow comparisons 
with theoretical predictions. The unfolding procedure corrects for effects collectively referred to as 'mass- 
bin migration’ (MBM) as described below. 

• Detector effects: the finite resolution of the detector, as well as lepton reconstruction efficiencies, 
deform the measured Z/y* line shape and the dependence of the asymmetry values on the dilepton 
mass with respect to what one would measure with an ideal apparatus covering the same kinematic 
range. 

• QED radiative corrections: radiative corrections [33], or real photon emission in the final-state 
(FSR), deform the shape of the dilepton invariant mass distribution. This deformation is partic¬ 
ularly pronounced below the Z peak. The events are moved from the Z peak (i.e. expected Afb 
positive and small) towards smaller values of invariant mass, significantly reducing the magnitude 
of the observed Afb in the region 66GeV< ma < mz . In the high-mass region (met > m z), the 
deformation due to radiative corrections is still present, but is reduced in magnitude. To account 
for these corrections, dileptons are unfolded to the pre-FSR state, referred to as ‘Born level’. 


The unfolding procedure is earned out using an iterative Bayesian unfolding method [34], as implemen¬ 
ted in the RooUnfold toolkit [35]. The response matrices are built from the PYTF1IA signal sample and 
the number of iterations (ten) is chosen in such a way as to optimize the result of closure tests on simu¬ 
lated samples. Additional checks are performed to ensure that the use of the PYTHIA LO generator for 
the unfolding does not bias the result. Since FSR is a significant correction, an alternative real-photon 
emission generator is investigated, using a simulated sample generated with SF1ERPA [36], which uses 
a module called PHOTONS++ [37] for higher-order QED corrections. The impact of NLO electroweak 
(EWK) corrections on the response matrix are estimated by reweighting the PYTHIA simulation to the 
prediction from the HORACE MC event generator and redoing the unfolding. In order to estimate NLO 
QCD effects on the A°g S values, a test is performed using a simulated POWHEG sample as pseudo-data 


2 Numerical values of all results are available in HepData [32], 


11 




« 1.2 

ctf 

e m 1 

Ea. 1 
< 0.8 
0.6 
0.4 
0.2 
0 

- 0.2 

-0.4 

- 0.6 

-2 


ATLAS 

is = 7TeV, 4.81b' 1 

CC electron 
p T > 25 GeV 
Ini < 2.47 






70 ^0 i 


=*= 


Data 

Xyfy/\ PYTHIA, Zff-tee 
A POWHEG, ZJf->et 


2x1Cr 


10 J 

m ee [GeV] 



(a) 


(b) 



Figure 3: Detector-level forward-backward asymmetry values as a function of the dilepton invariant mass for 

the (a) CC electron, (b) CF electron and (c) muon channels, after background subtraction. For the data, the black 
inner error bars represent the statistical component and the lighter outer error bars the total error (statistical and 
systematic added in quadrature). The boxed shaded regions for the MC expectations represent only the statistical 
uncertainty; theoretical uncertainties for MC are included in the systematic uncertainties on the data. The lower 
panel of each plot shows the pull value (A/cr) for each mass bin, where A is the difference between data and 
simulation and cr is the sum in quadrature of the data and simulation uncertainties. 
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Figure 4: Detector-level forward-backward asymmetry (A™g as ) values as a function of the dilepton invariant mass 
for the (a) CC electron, (b) CF electron and (c) muon channels in a narrow region around the Z pole, after back¬ 
ground subtraction. For the data, the black inner error bars represent the statistical component and the lighter outer 
error bars the total error (statistical and systematic added in quadrature). The boxed shaded regions for the MC ex¬ 
pectations represent only the statistical uncertainty; theoretical uncertainties for MC are included in the systematic 
uncertainties on the data. The lower panel of each plot shows the pull values (A/cr, as defined for figure 3). 
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and unfolding the asymmetry values using the PYTHIA-derived response matrices. These studies all 
show that any biases are much smaller than the present statistical uncertainties of the measurement. 

The systematic uncertainties on A°g S have contributions from the sources discussed in the following. 

• Unfolding uncertainty: estimated using a partially data-driven method. A set of weights is derived 
as a function of mu and cos ff cs to reweight the A j™ as values from simulation to the one observed in 
data. These weights are applied to the generator-level asymmetry values. The response matrix used 
in the unfolding is applied to the resulting values to fold and subsequently unfold them. Particular 
care is taken to make the matrices used for the folding and the unfolding statistically independent. 
The generator-level Afb dependence on mu is compared before and after the fold-unfold operation, 
and the difference is taken as an estimate of the uncertainty introduced by the unfolding. 

• Uncertainty due to finite size of the simulated event samples: the statistical uncertainty on the 
response matrices is propagated through the unfolding procedure. 

• Multi-jet background modelling: in the CC electron channel, the difference between the two back¬ 
ground estimation methods described in section 4.3 is taken as the systematic uncertainty and is 
found to be negligible with respect to other uncertainties. In the CF electron channel, this uncer¬ 
tainty is estimated by comparing templates based on different electron isolation requirements. For 
the muon channel, the impact of this background (and its uncertainty) is negligible. 

• Other experimental systematic uncertainties: these include the impact of pileup and of detector 
alignment, as well as energy/momentum scale and resolution, and trigger and reconstruction ef¬ 
ficiencies. The associated systematic uncertainties arc estimated following the prescriptions in 
Refs. [26, 28, 38]. The uncertainties related to energy scaling and resolution are among the largest 
contributions to the total error in all three channels, since they result in both a shifting and a broad¬ 
ening of the invariant mass peak, causing events to migrate between mass bins. 

• PDF uncertainties: the CT10 PDF set [39], which provides a reliable uncertainty estimate and is 
widely used in ATLAS, is also used here to estimate the PDF uncertainty. Its eigenvectors are 
used and the result quoted at 68% confidence level. For each error set, the MC signal sample is 
reweighted, the response matrices are recalculated and the unfolding is repeated. This contribution 
is found to be small when unfolding only mass-bin migration effects. 

The magnitudes of the systematic uncertainties on the A°g S values are summarized in table 1, for three 
invariant mass regions. Figure 5 shows the A°g S values obtained from leptons unfolded to Born level for 
all three channels. Expectations from PYTHIA and POWHEG are in good agreement with the measured 
values, as illustrated in the pull distribution at the bottom of each plot. 


5.1 Correcting for dilution 

A similar unfolding procedure is used to further correct the A°g S values to remove dilution effects, which 
occur when the wrong choice is made for the direction of the quark. The unfolding for dilution and 
the extrapolation from the detector acceptance to the full phase space are performed using the PYTHIA 
signal sample, where the description of the initial state allows a straightforward definition of the polar 
angle of the lepton with respect to the quark. The fully corrected asymmetry values for dileptons at 
the Born level Ap“' are shown in figure 6. The magnitude of the correction is larger than in the previous 
unfolding step. In addition, the contribution from the PDFs becomes the dominant systematic uncertainty. 
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Uncertainty 

CC electr 

66-70 GeV 

ons 

70-250GeV 

250-1000 GeV 

Unfolding 

~lxl0~ 2 

(2-5)xl0~ 3 

~4xl0 -4 

Energy scale/resolution 

~7xl0“ 3 

(0.5-2)xl0“ 3 

~2xl0“ 2 

MC statistics 

~5xl0“ 3 

(O.l-l)xlO -3 

(3-20)xl0 -3 

PDF 

~2xl0~ 3 

(l-8)xl0~ 4 

(0.7-3)xl0~ 3 

Other 

~lxl0~ 3 

(0.1-2)xl0 -3 

(5-9)xl0~ 3 

Uncertainty 

CF electr 

66-70 GeV 

ons 

70-250GeV 

250-1000 GeV 

Unfolding 

~2xl0~ 2 

(0.5-2)xl0“ 2 

- 

Energy scale/resolution 

~lxl0~ 2 

(0.5-7)xl0“ 2 

- 

MC statistics 

~lxl0~ 2 

(l-7)xl0~ 3 

- 

Background 

~3xl0" 2 

(0.5-l)xl0“ 2 

- 

PDF 

~4xl0“ 3 

(2-6)xl0~ 4 

- 

Other 

~lxl0“ 3 

(l-5)xl0~ 4 

- 

Uncertainty 

Muon) 

66-70 GeV 

> 

70-250GeV 

250-1000 GeV 

Unfolding 

~lxlO“ 2 

(l-4)xl0~ 3 

~5xl0~ 4 

Energy scale/resolution 

~8xl0“ 3 

(3-6)xl0~ 3 

~5xl0^ 3 

MC statistics 

~5xl0“ 3 

(O.l-l)xlO -3 

(2-30)xl0~ 3 

PDF 

~2xl0~ 3 

(l-8)xl0' 4 

(0.3-3)xl0^ 3 

Other 

~lxl0“ 3 

(0.5-l)xl0“ 3 

(3-10)xl0~ 3 


Table 1: Absolute systematic uncertainties on the A° 3s values, after unfolding for mass-bin migration. Approximate 
values in three invariant mass intervals are given. 


Good agreement is observed in general between the measured and predicted values. The muon channel 
measurement exhibits a discrepancy with respect to the PYTHIA prediction for masses above the Z boson 
mass, where the measured asymmetry is consistently larger than the prediction. This effect could not be 
explained by the analysis procedure and might be a feature of the simulation. 


6 Measurement of sin 2 6 lept 

e// 

The extraction of the effective weak mixing angle (sin 2 (r^ 1 ) from the detector-level asymmetry values 
(Aj^ as ) is presented here. 

Within the region of interest (0.218 < sin 2 6^ < 0.236) 17 MC simulated samples were generated with 
varying values of sin 2 The generator used for the templates is PYTHIA, which allows the value of 

sin 2 to be tuned without changing my. Within the range of the sin 2 9^ variations, the Z boson line 
shape remains unchanged in the generated samples. From the generator-level information in the samples, 
weights arc calculated, in bins of m// and cos6^ s , to transform the Aj™ as values to the ones expected 
(ApgO for a different value of sin 2 dj 1 ’ 1 . The reweighting technique is validated on simulated samples. 
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Figure 5: Forward-backward asymmetry (A°g S ) values as a function of the dilepton invariant mass for the (a) CC 
electron, (b) CF electron and (c) muon channels. Leptons are unfolded to Born level to account for mass bin 
migration, and the results are compared to truth-level MC information. For the data, the black inner error bars 
represent the statistical component and the lighter outer error bars the total error (statistical and systematic added in 
quadrature). The boxed shaded regions for the MC expectations represent only the statistical uncertainty; theoretical 
uncertainties are included in the systematic uncertainties on the data. The lower panel of each plot shows the pull 
values (A/cr, as defined for figure 3). 
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Figure 6: Forward-backward asymmetry (A“ r ) values as a function of the dilepton invariant mass for the (a) CC 
electron, (b) CF electron and (c) muon channels. Leptons are unfolded to Born level to account for mass bin 
migration and dilution effects are corrected. The measurement is extrapolated to the full phase space, and the 
results are compared to truth-level MC information. For the data, the black inner error bars represent the statistical 
component and the lighter outer error bars the total error (statistical and systematic added in quadrature). The 
boxed shaded regions for the MC expectations represent only the statistical uncertainty; theoretical uncertainties 
are included in the data systematic uncertainties. The lower panel of each plot shows the pull values (A/cr, as 
defined for figure 3). 
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For each channel, the A^ values obtained from the reweighted datasets are compared to those obtained 
from the data, using a x 2 test over the mass range 70-250 GeV, taking statistical and systematic uncer¬ 
tainties into account. The mass range has been optimized for the maximum sensitivity and stability of 
the measurement. A parabola is fitted to the resulting distributions of x 2 values for each channel inde¬ 
pendently. The minimum of the parabola yields the sin 2 result. The x 2 /ndf value at the minimum 
is 22.4/16 for the CC electron channel, 21.9/16 for the CF electron channel and 22.6/16 for the muon 
channel. The fit results are found to be stable with respect to the invariant mass range over which the 
template comparisons are performed, as well as with respect to the sin 2 9^ range over which the x 2 is 
minimized. 

As discussed in section 5, the use of a LO generator and a specific implementation of the real photon 
emission in the final state does not bias the unfolded A J™ as values. In order to assess the impact of these 
potential sources of systematic effects on the templates used to extract the sin 2 9^ value, additional tests 
are performed. Effects related to the modelling of the real photon emission are tested with SF1ERPA, and 
are found to be negligible. The effects of NLO QCD corrections are investigated further in the context 
of the sin 2 d^j measurement by comparing LO and NLO predictions of the Afb vs. mass distributions 

calculated with MCFM [18]. Differences are propagated through the extraction of sin 2 6^ and the res¬ 
ulting variation of the weak mixing angle is treated as an additional systematic error. The effects of NLO 
EWK corrections are found to be small compared to the rest of the uncertainties, but are accounted for as 
an additional systematic uncertainty on the final result. 

The systematic uncertainties already described for the A™g as values are also estimated for the sin 2 8^ 
measurement. As the background is small in all channels of the sin 2 6^ extraction, a simple, but slightly 

conservative, approach is used to obtain its uncertainty. The sin 2 measurement is repeated without 
the subtraction of the background, and the result is compared to the baseline measurement, which has 
the background subtracted. The uncertainty on sin 2 0^ 1 from the background is taken to be 10% of 
the observed difference, to take into account the uncertainties on the cross sections of the background 
components known with least precision [40] . 


6.1 Impact of PDFs on the sin 2 measurement 

This measurement is sensitive to the PDFs describing the flavour composition of the initial state, since 
the Afb values depend on the flavour and charge of the initial partons. In addition, the u v and d v valence 
quark distributions have an impact on the measurement due to differences in the weak couplings, while 
dilution effects introduce a dependence on the sea-to-valence quark ratio within the accessible Bjorken-v 
range. 

The ATLAS measurements of inclusive W and Z boson production [41] are sensitive to the same effects 
and indicate that some of the existing PDF sets may not provide a good description of the data. 

For this measurement of sin 2 a LO version of the ATLAS-epWZ NNLO and NLO fits was prepared. 
In the following, the PDFs extracted from these fits are called ATLAS-epWZ12 LO PDFs. These fits 
include inclusive W + , W~ and Z production data at \[s = 1 TeV [41], together with the combined F1ERA 
data on inclusive neutral- and charged-current interactions from e + p and e~p scattering [42]. The set¬ 
tings for the fit are kept the same as those for the NLO and NNLO fits as much as possible. The main 
differences between the LO fit and the higher-order fits are that the gluon PDF distribution has a simple 
parameterization which requires it to be positive definite at all scales, and that the value of the strong 
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coupling constant a s (mz) is set to be higher than the one for the NLO or NNLO fits. A value of a s (mz) 
of 0.130 is chosen with a scanning procedure, such that it yields the best level of agreement (x 2 ) between 
data and the fit result. This value is consistent with that used by other LO PDF sets. These PDFs are 
available in LHgrid formats [43]. 

6.2 Results for sin 2 6 lcpt 

e// 

The measured values of the weak mixing angle obtained with ATLAS-epWZ12 LO PDF are shown in 
table 2. Of the three channels, the CF electron channel has the lowest statistical uncertainty, despite 
having the fewest selected events, as detailed in section 4.3. This is because the Apg as values in this 
channel are less affected by dilution due to the larger average rapidity of the dilepton system compared to 
the other two channels. Details of the main sources of systematic uncertainty on the result in each channel 
are given in table 3. 

The sin 2 9 ll i p , 1 measurements from all three channels are combined. Given the total covariance mat- 

e// 

rix C, the weights assigned to the three measurements for the combination are calculated using uy = 
Yik Q 1 / Yj jk C~ k , with the indices i, j and k representing the three measurements (CC electron, CF elec¬ 
tron, muon). The total error is given by cr 2 = (W T C _1 ) _1 VP, with W T = (1,1,1). Uncertainties due to 
energy/momentum scale and resolution are treated as completely uncorrelated across channels. Since the 
energy scale and resolution for the CF electron channel is dominated by the forward calorimeter perform¬ 
ance, this is a good approximation for the correlation of the CC and CF electron channels. The systematic 
uncertainties due to the MC statistical uncertainty are also treated as fully uncorrelated. Theory-related 
systematic uncertainties, such as those associated with PDFs and higher-order EWK and QCD correc¬ 
tions, are treated as fully correlated across channels. All the remaining uncertainties are treated as fully 
correlated across the channels to which they are applicable. The central value of the combination is found 
to be stable when varying the magnitudes of the main correlated uncertainties. Following this procedure, 
the combined (CC+CF) electron result, as well as the electron-muon combination, is shown in table 2. 
The systematic uncertainty on the combined result is dominated by the PDF uncertainty (±0.0009), which 
is calculated using the variations provided in the ATLAS-epWZ12 LO PDF eigenvalue set, similar to the 
calculation described in section 5. The contributions to the total systematic uncertainties on the results 
are included in table 3. 

The results obtained with the ATLAS-epWZ12 LO PDF set have been compared to those obtained using 
other leading order PDF sets. The PDF sets ATLAS-epWZ12 LO and F1ERAPDF1.5LO [44] yield very 
similar results, the ATLAS-epWZ12 LO result being larger by 1 x 10~ 4 (with a PDF uncertainty of 
9x 10 -4 ). The MSTW2008LO set produces a downward shift in the resulting sin 2 which is significant 

(A sin 2 6 e ^ pt = -2 x 10 3 ). However, MSTW sets are known to give a poor description of the ATLAS W 
and Z data [41]. The CT10 PDF set (which is NLO) was also tested, and yields a sin 2 6^ pt value which is 
smaller than, but compatible with, the results obtained with the ATLAS-epWZ12 LO PDF (A sin 2 6^ pt = 
-8 x 1(T 4 ). The PDF analysis of the ATLAS data also suggests an increased strange-quark sea density 
( s/d ~ 1) compared to other PDF sets [15]. In order to probe the sensitivity of this measurement to 
the enhanced strange-quark density, a dedicated PDF set was prepared with a suppressed strange sea 
corresponding to uncertainty, which indicates a low sensitivity of the measurement to this effect. This is 
due to the fact that the sea composition only affects the measurement through the dilution, which, to first 
approximation, does not change the position of the minimum of they' 2 in the template fits. 


19 


As a cross-check of the measured weak mixing angle and the unfolding procedure, sin 2 9 ^ i s extracted 
from the unfolded particle-level asymmetries. Similarly to the procedure described in the beginning of 
section 6 , samples with various values of sin 2 8^, generated using MCFM v 6.8 [18] at LO in perturbative 
QCD and interfaced to APPLGRID v 1.4.69 [45] using the ATLAS-epWZ12 LO PDF set, are compared 
to the measured particle-level asymmetries using a y~ test. The correlations between the systematic 
uncertainties are included in the y 2 calculation. A parabolic fit to the y 2 distribution yields a value of 
sin 2 6 c jj in good agreement with the detector-level extraction. 

There are several other measurements of sin 2 6^ at the m? scale. The results from this analysis, as well 
as the results from the other collider experiments, are summarized in table 4 and shown in figure 7. These 
include measurements from LEP, SLC, the Tevatron, and the LF1C, as well as the results of the PDG 
global fit [46]. Differences with respect to the combined LEP/SLC measurements are also displayed. The 
combined result of this analysis agrees within 0.1 cr with the most precise leptonic asymmetry measure¬ 
ment, and within 1.2 cr of the sin 2 8 le ?} value extracted from the LEP A^:!’ measurement. The combined 
result is in agreement with the PDG global fit at the level of 0.6 cr. 



sin 2 fl lept 

Sm °eff 

CC electron 

CF electron 

Muon 

El. combined 

0.2302 ± 0.0009(stat.) ± 0.0008(syst.) ± 0.0010(PDF) = 0.2302 ± 0.0016 

0.2312 ± 0.0007(stat.) ± 0.0008(syst.) ± 0.0010(PDF) = 0.2312 ± 0.0014 

0.2307 ± 0.0009(stat.) ± 0.0008(syst.) ± 0.0009(PDF) = 0.2307 ± 0.0015 

0.2308 ± 0.0006(stat.) ± 0.0007(syst.) ± 0.0010(PDF) = 0.2308 ± 0.0013 

Combined 

0.2308 + 0.0005(stat.) ± 0.0006(syst.) ± 0.0009(PDF) = 0.2308 + 0.0012 


Table 2: The sin 2 measurement results in each of the three studied channels: electron central-central, electron 
central-forward and muon. Results of the statistical combination of both electron channels and all three channels 
are shown as well. 


Uncertainty source 

CC electrons 
[ 10 ~ 4 ] 

CF electrons 
[ 10 - 4 ] 

Muons 

[ 10 - 4 ] 

Combined 

[ 10 - 4 ] 

PDF 

10 

10 

9 

9 

MC statistics 

5 

2 

5 

2 

Electron energy scale 

4 

6 

- 

3 

Electron energy resolution 

4 

5 

- 

2 

Muon energy scale 

- 

- 

5 

2 

Higher-order corrections 

3 

1 

3 

2 

Other sources 

1 

1 

2 

2 


Table 3: Contributions to the systematic uncertainties on the sin 2 8^ values extracted from the three analysis chan¬ 
nels and on the combined result. Null entries (denoted by correspond to uncertainties that do not apply to a 
specific channel. Higher-order corrections include NLO QCD and NLO EWK contributions. Other sources include 
the effect of pileup, background uncertainties, lepton trigger/reconstruction/identification efficiency uncertainties, 
muon momentum resolution and effects of detector misalignment. 
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A/cr 

A/cr 


sin 2 9 

e// 

(wrt LEP+SLC) 

(wrt ATLAS) 

CMS [6] 

0.2287 ± 0.0032 

-0.9 

-0.6 

DO [5] 

0.23146 + 0.00047 

-0.1 

0.5 

CDF [4] 

0.2315 + 0.0010 

-0.03 

0.4 

LEP, [3] 

0.23221 +0.00029 

- 

1.2 

LEP, [3] 

0.23099 ± 0.00053 

- 

-0.1 

SLC, A lr [3] 

0.23098 ± 0.00026 

- 

-0.1 

LEP+SLC [3] 

0.23153 + 0.00016 

- 

0.6 

PDG global fit [46] 

0.23146 + 0.00012 

-0.4 

0.6 

ATLAS 

0.2308 ± 0.0012 

-0.6 

- 


Table 4: Comparison of the results of this analysis with other published results for sin 2 0^. The comparison 
includes the most precise measurements from LEP and SLC, and the results from the leptonic sin 2 meas¬ 
urements from the hadron collider experiments CMS, DO, and CDF. Also shown are the values of sin 2 9^ from 

the LEP+SLC global combination (which includes all sin 2 9^ measurements performed at the two colliders) and 
from the PDG global fit. Each A/cr column shows the difference between the result and the quoted reference value, 
divided by the quadratic sum of the associated uncertainties. 


6.3 Determination of 

As described in section 1, the forward-backward asymmetry around the Z pole for the muon channel, 
Apg, can be expressed in terms of the muon and quark asymmetry parameters A fJ and A q , c.f. eq. 4 and 
ref. [3]. Both are functions of the vector and axial-vector couplings of the quark/muon, 


Aq/fj - 


2gfgf 

(5>f ) 2 + (af ) 2 


2 gflgf 

l + < dthf)-' 


(5) 


The asymmetry parameters are related to the flavour-dependent weak mixing angle sin 2 6^ by 


gflgT = i - 4|G, /P I sin 2 e $ f , (6) 

where Q q /^ is the quark (q) or muon (ju) charge. The measurement of the forward-backward asymmetry 
can be interpreted as a determination of when assuming sin 2 6^ = sin 2 9^ = sin 2 which is 
valid within an uncertainty of 1.5 x 10 4 [3]. Given the value of sin 2 9^, A fl is small compared to A q and 
the uncertainty on sin 2 9^ (~ 1.5 x 10 -3 ) plays a more important role in eq. 4 than the error introduced 
by the assumption about sin 2 0^ mentioned above. With the present precision, it is therefore possible to 
use the value of sin 2 P' obtained from the muon asymmetry measurement to derive 


2(1 -4sinV;P l ) 

1+(1-4 sin 2 0^)2' 
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Figure 7: Comparison of the results of this analysis with other published results for sin 2 0^ p '. This includes the 
most precise measurements from LEP and SLC, and the leptonic sin 2 8^ measurements from the hadron collider 
experiments CMS [6], DO [5], and CDF [4], Also shown are the values of sin 2 6^ from the LEP+SLC global 

combination [3] (which includes all sin 2 O^jy measurements performed at the two colliders) and from the PDG 
global fit [46]. The vertical dotted line shows the central value of the ATLAS combined measurement reported 
here, while the vertical dashed line represents that of the current PDG global fit [46]. 

The uncertainties on A M are propagated from the uncertainties on sin 2 O^ff- result is 

= 0.153 ± 0.007(stat.) ± 0.006(syst.) ± 0.007(PDF) = 0.153 ± 0.012(tot.), 

which is of similar precision and in good agreement with the measurement from e + e~ collisions of 0.142+ 
0.015 [3]. It is worth stressing, however, that the determination of A t , in the LEP/SLD results is based 
entirely on asymmetry measurements in the different lepton final states without any assumptions on other 
Af, whereas the determination of A p presented here uses the Standard Model prediction of A q . 


7 Conclusions 

The forward-backward asymmetry in electron and muon pairs from Z/y* decays is measured using the 
7 TeV pp LHC collision data recorded with the ATLAS detector in 2011 corresponding to an integrated 
luminosity of 4.8 ftr 1 . The data are analysed over a range of dilepton invariant masses from 66 GeV to 
1000 GeV in the central-central electron and muon channels, and up to 250 GeV in the central-forward 
electron channel. The latter includes events where one electron is reconstructed in the forward pseu¬ 
dorapidity range (2.5 < \r/\ < 4.9). 
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The forward-hack ward asymmetry is measured separately for the three channels as a function of the 
dilepton invariant mass and unfolded for detector effects and final-state radiation. Additionally, a leading- 
order interpretation which accounts for the effects of dilution and full detector acceptance is presented. 
The resulting Afb values arc found to be in agreement with the corresponding Standard Model predic¬ 
tions. 

The detector level asymmetry values arc used to extract the value of the leptonic effective weak mixing 
angle, sin 2 separately for the three data samples using ax 2 minimization method. The results are in 
good agreement with each other and with measurements at e+e~ colliders, at the Tevatron and by CMS at 
the LHC. 

Results from the electron and muon final states arc combined, yielding 

sin 2 ej 1 * = 0.2308 + 0.0005(stat.) ± 0.0006(syst.) ± 0.0009(PDF) = 0.2308 ± 0.0012(tot.). 

The dominant uncertainty comes from knowledge of the PDFs. 

The result from the muon channel, when converted to the asymmetry parameter d /( , yields 

= 0.153 ± 0.007(stat.) ± 0.006(syst.) ± 0.007(PDF) = 0.153 ± 0.012(tot.), 
which is in good agreement with the best previous measurements. 
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